Green's function for magnetically incoherent interacting electrons in one dimension 
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Using a path integral approach and bosonization, we calculate the low energy asymptotics of the 
one particle Green's function for a "magnetically incoherent" one dimensional strongly interacting 
electron gas at temperatures much greater than the typical exchange energy but much lower than the 
Fermi energy. The Green's function exhibits features reminiscent of spin-charge separation, with 
exponential spatial decay and scaling behavior with interaction dependent anomalous exponents 
inconsistent with any unitary conformal field theory. We compute the tunneling density of states 
at low energies and find that it is a power law in energy with exponent l/{Ag) — 1, where g is the 
Luttinger interaction parameter in the charge sector. The underlying physics is made transparent 
by the simplicity of the approach. Our results generalize those of Cheianov and Zvonarev [Phys. 
Rev. Lett. 92, 176401 (2004)]. 
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Interacting one dimensional electron systems have 
proven to be especially rich in their physics, most no- 
tably because of the universal low energy properties that 
are present when the interactions are not too strong, 
the so-called Luttinger Liquid (LL) state whose exis- 
tence in nature is now establishedQ- It is well knownQ 
that the LL state exhibits decoupled spin and charge de- 
grees of freedom (spin-charge separation) with distinct 
spin and charge velocities for the respective collective 
modes. Characteristic of the LL state are power-law de- 
cays of various correlation functions, notably the single- 
particle Green's function, which is suppressed relative to 
the naive expectation of Fermi liquid theory. As one con- 
sequence, the local tunneling density of states vanishes 
in a power-law fashion as the chemical potential is ap- 
proached. All low-energy properties of the LL state can 
be understood from bosonization, which relates them to 
correlators in a simple free boson unitary conformal field 
theory (CFT). 

In the regime of strongly interacting, very low den- 
sity electrons, different physics is to be expected when 
the assumptions of LL theory break down. The distinc- 
tion between low and high densities is often quantified 
by the parameter = (nas)"^ where n is the average 
electron density and as is the Bohr radius specific to the 
material. When rg ^ 1 the spacing between electrons is 
large compared to ob and the potential energy dominates 
the kinetic energy, driving the system towards a Wigner 
crystal. When such strong interactions are present, it 
becomes difficult for electrons to exchange their position 
since they must effectively tunnel through one another, 
leading to an exponentially small0, Q (in r^) exchange 
energy, J. It then becomes quite easy to reach the mag- 
netically incoherent regime where the exchange energy 
is much less than the temperature: J <^ T. In this 
Letter we compute the low energy asymptotics of the 
1-d Green's function for general interactions in the limit 



Ts 3> 1 when the temperature is still much less than the 
Fermi energy: J <C T ^ Ep- We find that the Green's 
functions exhibits exponential decay in the spin sector 
and power law decay in the charge sector characterized 
by interaction dependent anomalous exponents (which 
do not correspond to any unitary conformal field the- 
ory). Our results generalize the results of Cheianov and 
Zvonarev[Q| (CZ) and the zero field results of Berkovich0 
and are obtained in a much simpler and more physically 
transparent manner. 

Due to breakthroughs in materials technology that 
allow unprecedented exploration of clean 1-d quantum 
wires "s^ theoretical interest in such systems has also been 
renewed 9J . The regime J <^T ^ Ep was recently con- 
sidered by Matveevpj where he showed there is a drop in 
the conductance of a one channel wire from 2e^//i to /h 
when J drops below T. In the same regime CZ^] have 
computed the low energy asymptotics of the one particle 
Green's function assuming infinitely strong zero range in- 
teractions between electrons. Their tour-de-force Bethe- 
ansatz based calculation, however, does not provide clear 
insight into the physics. Moreover, their results are spe- 
cific to the special features of the Hamiltonian they con- 
sidered. Our calculation of the 1-d Green's function in 
the regime J <^T <^ Ep is completely general and sim- 
ple enough to highlight the physical origins of the non-LL 
features. 

We note that the physics of the J <^T <^ Ep regime 
is of broad and general interest as it may play some role 
in other systems not fully understood, such as the two di- 
mensional metal-insulator transition and the 0.7 anomaly 
in quantum point contacts. As such, it is a worthwhile 
endeavor to elucidate the basic physics of this regime. 

MODEL: We assume that the electrons are confined 
to one dimension and experience predominantly repul- 
sive and spin independent interactions, i.e. ones which 
can be written solely in terms of the local electron den- 
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sity n{x) = J2a=-\ i ^ct(^)V'o-(^)- Here ipaix) is the field 
operator for an electron at position x with spin a. In the 
discussion that follows, we will not require an explicit 
Hamiltonian, but rather make use of an eff'ective low- 
energy theory that contains renormalized parameters, Vc 
and g, depending on the microscopic interactions in an 
unspecified way. Despite our interest in the magnetically 
incoherent regime J <C T, we will see that the essen- 
tial parameters describing the effect of interactions may 
be obtained from the charge sector action of LL theory 
which obtains at the lowest temperatures T <^ J: 
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where Vc is the velocity of charge excitations, g is the pa- 
rameter of the low energy theory measuring the strength 
of electron interactions and the charge fields 9 and </> are 
the fields appearing in the low-temperature bosonized 
version of the electron operator. The charge fields are 
defined via 9^ +61^9 and ((/)| + (?!)x)/2 = 0. The 9 field 
is related to the particle density fluctuations through the 
familar relationship n{x) — ^dx9[x) and e""^^^^ annihi- 
lates/creates a particle at x. We choose to do all our 
calculations in imaginary time, r = it. The Green's func- 
tions of interest (retarded, advanced, etc.) can then be 
computed by the appropriate analytical continuation to 
real time. 

RESULTS: Wc compute the single particle Green's 
function 



g{x,T) = {^k^{x,T)i>\{{),0)) 



(2) 



in the limit J <^ T Ep hy first averaging over spin 
configurations and then computing the low energy charge 
dynamics using Due to the spin rotational invari- 
ance of the electron interactions, the Green's function 
for spin down electrons coincides with We find (for 
a; > 0, r > ^ 00), 



VnT — IX 
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VcT + ix 



where C is an undetermined constant (CZ determine it 
for the special case of infinite strength zero range interac- 
tions in Ref. |^). Here we follow CZ and define a spinless 
Fermi wavevector kp = T^n, where n is the average den- 
sity of electrons. The phases (p^ are given by 
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(4) 

The power law decay of Q is characterized by the 
anomalous exponent 




(5) 



It is clear that the Green's function Q does not fit 
into the usual LL paradigm of correlation functions with 
power law decay because of the exponential decay with 
distance. We will show that this exponential decay is 
simply a result of spin averaging when J <C T. The ex- 
ponents characterizing the power law decay of the charge 
sector are anomalous because for certain values of 5, 
g = 1/2 for example, A^^i = -\ {^f < 0, which, if 
interpreted as arising from CFT, implies non-unitarity in 
the charge sector 6]. We will show the anomalous power 
law decay of © comes from density fluctuations in the 
charge sector after averaging over the spin degrees of free- 
dom. 

We also compute the x = Green's hmction 
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from which the low frequency spectral function can be 
computed 
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All our results recover the results of CZ in the spe- 
cial case they considered of infinite repulsive local {6- 
function) interactions, which is just the infinite U limit 
of the Hubbard model. The CZ limit corresponds [loj to 
g ~ 1/2 from which all their results can be recovered by 
plugging this value of g into Eqs. |0J-|2I).(CZ ignore the 
y^\n{w) factor in our results agree exactly at g=l/2.) 

DERIVATION OF RESULTS: We study the single- 
particle Green's function for a general (finite range and 
strength of interactions), strongly interacting 1-d elec- 
tron system at finite temperature: 
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Here Z = Tt [e~^^] is the partition function and /3 is the 
inverse temperature. Our results are based on the first 
quantized path integral representation of Gix, r) and Z 
in imaginary time, r = it. As is well-known, Z is ob- 
tained as an integral over world lines (paths) of up and 
down spin electrons, with periodic boundary conditions 
in imaginary time allowing for permutations of identical 
particles, i.e. the final positions of electrons of a given 
spin polarization must be a permutation of their initial 
positions (up and down electrons may not be permuted 
as they are not identical). Each such configuration is 
weighted by (— l)^T+fie~'^=, where Se is the Euclidean 
action describing the "deformation" of the world lines 
and (—I)-'''' is the sign of the permutation of spin a elec- 
trons. The numerator of Eq. ((SJ is a similar integral over 
paths, with a spin-up electron world line inserted at (0, 0) 
and another removed at (x,t). To compute the sign of 
the permutation, one should treat the path terminating 
at {x,t) as continuing from (0,0). 



A crucial element of our analysis is a non-crossing con- 
dition. Due to the Pauli principle, world lines of the same 
spin electrons can be treated as non-crossing irrespective 
of interactions. Moreover, the condition J ^ T precisely 
corresponds to the absence of crossings (exchanges) of 
opposite spin world lines {1/J:s>(3 gives the typical dis- 
tance in imaginary time between exchange events) . Phys- 
ically, this is due to the large Coulomb repulsion of elec- 
trons. Thus the topology of the paths is identical to that 
of spinless fermions or infinitely repulsive spinless bosons. 
Some world line trajectories are shown in Fig. ^ 

We choose to compute the trace in Eq. (jHl by first 
summing over all possible spin assignments for each set 
of world lines. This gives an effective weight for the re- 
maining sum over the world line configurations. This 
weight includes two contributions: (1) the Euclidean ac- 
tion factor e~^^, and (2) a magnetic/statistical factor ob- 
tained from the spin sum and permutation signs, which 
is discussed below. Given the non-crossing topology and 
spin-independent interactions, the Euclidean action is ex- 
pected to be well-approximated (for T <^ Ep) by the 
low-energy effective form appropriate to spinless (in gen- 
eral interacting) fermions/bosons, which justifies Eq. l|T]). 
Moreover, for J <C Ep, we expect this "charge sector" 
action to be identical to the ultimate low-energy charge 
action of LL theory valid for T <C J. 

Next consider the spin sum. The limit J <^T implies 
that for a fixed set of world lines all spin configurations 
are computed with equal weight. In order for a particle 
created at spacetime coordinate (0,0) to be annihilated 
at (x, r) a large distance away, the world line must "wrap 
around" the imaginary time inverval to /3 a large num- 
ber of times. From Fig. ^ it is clear that to stay in the 
low energy sector, all intermediate world lines must have 
the same spin. Any configuration of world lines without 
all the intermediate lines having the same spin would vi- 
olate the non-crossing condition. When the spin average 
is taken in the trace in ||HJ), only the term with all inter- 
mediate worldlines having spin up will contribute to the 
low energy Green's function. For N intermediate world 
lines, the probability that all iV -I- 1 spins are the same 
is 2^^. After an electron added at (0,0) reaches (a;,r), 
it has been permuted through iV electrons of the same 
spin and therefore the Green's function picks up a factor 
(-1)^. 

After spin averaging for J ^ T, the Green's function 
therefore becomes 

g{x,T) - (2-^(^'^)(-l)^(^'^)e^('^(^'^)-'^(0'0))) , (9) 

where now the average is taken at T ^ in the charge 
sector. The term e*('^(^'^)-'^(o>o)) creates a world line at 
(0,0) and annihilates it at (cc, t). (Note that at this stage 
in the calculation the dynamics after spin averaging be- 
come effectively spinless as discussed earlier.) All the 
effects of statistics and spin are encapsulated in the first 
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Particle world lines r »1 




X X 




X X 

FIG. 1: World lines for a strongly interacting 1-d electron 
gas at J « T. Particle trajectories in space and imaginary 
time are shown as curved lines. The dashed lines represent 
the world line paths for creating a particle and removing it 
for large x = xj — Xi,T = tj — Ti, ((a) and (b)), and for 
X = Q,T ((c) and (d)). The solid lines represent trajectories 
of other particles. Due to the large action cost associated 
with the trajectories in (b) , at low energies a process like that 
shown in (a) where world lines "wrap around" from r = /3 
to r = will dominate. Such a process, however, requires 
that all dashed world lines have the same spin. For J <^ T 
this occurs with probability as discussed in the text. 

Fig.(c) shows a contribution to the k — 1 term of Eg. 11131 . 
however processes like that shown in (d) dominate and yield 
the Green's function 

two terms inside the average. The number of electrons is 
related to the 9 field via 

N{x, t) = fix -f i {0{x, t) - 61(0, 0)) , (10) 

TT 

which expressed the number of electrons in a distance x in 
terms of the average density and a small fluctuating piece 
expressed in terms of the 9 fields. Using Eq. H10|) . and 
writing (—1)^ = Re[e"^], the simplest form correct for 
integer N (the harmonic approximation violates this) and 
consistent with the requirement that Q(x^t) be real and 
even in x. The Green's function can now be expressed as 

g = g+ + g^, with g^ = [g+]* and 

g+{x,T) ^ ^-k.'.ci-^^k,.^^-^-^(e(.,r)-e(m) 



4 



where the first two terms come from the exponentiation 
of the average density and we have used kp = i^n. This 
clearly identifies the exponential decay of the first term 
as coming from spin averaging and the oscillatory second 
term as coming from Fermi statistics. 

We now compute the part of the Green's function com- 
ing from fluctuations in the charge sector using the action 
(P). Making the definitions, <i>(x,T) = 4>(,x,t) - (/)(0,0) 
and 6(a;, r) = 6[x, r) — 6'(0, 0), we use the Gaussian ac- 
tion to move the averages to the exponent, 

xe-i{*^)e-(i+'^)<*e) _ (12) 

Standard computations from Eq. give (O^) — 
g\n{x^ +vIt^), ($2) = ^ln(a;2 + vIt"^) and ($8) = 

1 (^i^TTfi) ■ Substituting these values into Eq. (O, 
the anomalous exponent given in Eq.ljSJ is obtained. An 
additional phase factor comes from the — — ((0^) + 

(*e)) = -^^(g\n{x^+vlT^) + \\n[^)) piece 
in the exponent. Combining this with the complex con- 
jugate yields Eq. ©. 

The preceding calculation brings out the essential 
physics of the Green's function Q : The spin averaging is 
responsible for the exponential decay of the Green's func- 
tion and imposes a constraint on the world line configura- 
tions that contribute to it. Fermi statistics is responsible 
for the oscillatory terms. Treating Gaussian fluctuations 
about this constraint in the charge sector results in the 
power law decay with generalized anomalous exponents 
depending on the interaction parameter g. 

Having discussed the spatial asymptotics of ^(a;, r), for 
X, T ^ oo, we now turn our attention to Q{Q, r) which will 
allow us to compute the low energy tunneling density of 
states at a point. Unlike the situation with a; oo, 
computing the Green's function at x = forces us to 
consider the discreteness in the number of world lines 
that may "bend' in between (0,0) and (0,t): 

oo 

g{Q,T) - ^ 2-l'^-|(-l)'=((5(]V(0,r) -fc)e'('^(°'^)-'*("^o») 

k— — oo 

1 °° 2 .2 

^ ^ 2-l'=l(-l)'=e-3^^fe , (13) 

where the result is recovered by noting that the sum 
depends only weakly on r and ranges between 2/3 and 
1. Fourier transforming H13() into frequency space imme- 
diately gives (|7|). 

DISCUSSION: We have computed the low energy 
asymptotics of the one particle Green's function g{x, t) 
in the limit of » 1 and J <^ T <^ Ep for arbi- 
trary interactions. We find the correlation function does 



not fit the usual LL form. Instead, the spin averaging 
present in the Green's function for J <C T results in an 
exponential decay of the Green's function with distance. 
The low energy behavior of the charge sector still shows 
a power law decay due to Gaussian fluctuations, but 
with interaction dependent anomalous exponents. The 
low frequency tunneling density of states (proportional 
to the spectral function) also shows interesting behav- 
ior depending on g: A{ijo) oc uj'^/'^^s)-'^ ^ which shows a 
crossover from a power law divergence for g > 1/4, to a 
pseudo-gap when g < 1/4. The divergence at low ener- 
gies can be understood as coming from the infinite spin 
degeneracy when J <^ T . As the interactions increase 
{g decreases) the suppression of the density of states in 
the charge sector overwhelms the spin degeneracy to re- 
cover the power law suppression familiar in LL theory. 
Finally, we note that Fourier transforming the Green's 
function (at low frequency) into momentum space will re- 
sult in broad peaks of width ~ kp centered at k « ±fci?. 
Note that the "Fermi momentum" hp appearing here dif- 
fers by a factor of 2 from the usual Fermi momentum: 
kp — 7rn/2 = kp/2. Hence, as one moves from Ts « 
to rs 3> 1 and the regime J T <^ Ep is reached, 
one expects to see delta function-like peaks present at 
small rs to broaden to width kp and the centers to 
shift from ±7rn/2 to ±7rn creating a broad double- lobed 
structure in momentum resolved tunneling when 3> 1 
and J <^ T <^ Ep. We hope this work will inspire new 
ideas in other systems where the physics discussed here 
may play a role. 
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